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In this paper, we consider the bosonic t- J model, which describes two-component hard-core bosons 
with a nearest-neighbor (NN) pseudo-spin interaction and a NN hopping. To study phase diagram 
of this model, we derive effective field theories for low-energy excitations. In order to represent 
the hard-core nature of bosons, we employ a slave-particle representation. In the path-integral 
quantization, we first integrate our the radial degrees of freedom of each boson field and obtain 
the low-energy effective field theory of phase degrees of freedom of each boson field and an easy- 
plane pseudo-spin. Coherent condensates of the phases describe, e.g., a "magnetic order" of the 
pseudo-spin, superfluidity of hard-core bosons, etc. This effective field theory is a kind of extended 
quantum XY model, and its phase diagram can be investigated precisely by means of the Monte- 
Carlo simulations. We then apply a kind of Hubbard- Stratonovich transformation to the quantum 
XY model and obtain the second-version of the effective field theory, which is composed of fields 
describing the pseudo-spin degrees of freedom and boson fields of the original two-component hard- 
core bosons. As application of the effective- field theory approach, we consider the bosonic t-J model 
on the square lattice and also on the triangular lattice, and compare the obtained phase diagrams 
with the results of the numerical studies. We also study low-energy excitations rather in detail in 
the effective field theory. Finally we consider the bosonic t-J model on a stacked triangular lattice 
and obtain its phase diagram. We compare the obtained phase diagram with that of the effective 
field theory to find close resemblance. 

PACS numbers: 67.85.Hj, 75.10.-b, 03.75.Nt 



I. INTRODUCTION 

In recent years, systems of cold atoms have attracted 
interest of many condensed- matter physicists. In par- 
ticular, cold atoms put on an optical lattice sometimes 
regarded as a "simulator" to study canonical models 
of strongly correlated electron system^. These sys- 
tems are quite controllable and impurity effects are sup- 
pressed. From this point of view, we have started study 
of the bosonic t-J model^^^, which is a bosonic counter- 
part of the t-J model for the high- Tc superconducting 
phenomena"^. The bosonic t-J model (B-t-J model) de- 
scribes two-component bosonic cold atoms in an optical 
lattice with strong inter and intra-species repulsions that 
can be controlled in the experiments. 

Our previous studies on the B-t-J model mostly em- 
ployed numerical Monte-Carlo (MC) simulations to clar- 
ify its phase diagram at finite temperature (T)^. In con- 
trast to the fermionic t-J model, the numerical study can 
be done without any difficulties for some cases of the B-t- 
J model. In this paper, we shall derive effective fields the- 
ories of the B-t-J model and study them by both analyti- 
cal and numerical methods. We mostly focus on the B-t-J 
model on square and triangular lattices at T = and the 
phase diagram for quantum phase transitions. The ob- 
tained results are compared with the previous findings 
and numerical study of finite-T systems. 

This paper is organized as follows. In Sec. II, we shall 
introduce the B-t-J model and derive first- version of an 
effective field theory by using the path-integral meth- 
ods. The local constraint is faithfully treated by using 



the slave-particle representation. By integrating out am- 
plitude modes of the slave particles, we obtained an ex- 
tended quantum XY model, which describe Bose conden- 
sation of two atoms and pseudo-spin degrees of freedom. 
Numerical study by the Monte-Carlo (MC) methods can 
be performed without any difficulties and the phase dia- 
gram is obtained for the B-t-J model on both the square 
and triangular lattices in Sec. III. In Sec. IV, we shall de- 
rive second- version of the effective field theory from the 
extended quantum XY model by introducing collective 
fields for the two bosons and the pseudo-spin degrees of 
freedom. By using this effective field theory, we obtain 
phase diagram of the B-t-J model and verify the consis- 
tency of the results. Furthermore we study low-energy 
excitations, i.e., the Nambu-Goldstone bosons in various 
phases in the phase diagram. In Sec.V, we numerically 
shall the B-t-J model in the stacked triangular lattice at 
finite T. Obtained phase diagram has a similar structure 
of the model in the triangular lattice at T = 0. Section 
VI is devoted for conclusion and discussion. 

II. BOSONIC t-J MODEL AND DERIVATION 
OF EFFECTIVE MODEL 

Hamiltonian of the B-t-J model, which will be studied 
in this paper, is given as^^, 

H,, = -Y, + h\h, + h.c.) + J, ^ 5f 5| 

{i,3) {ij) 

+jJ2iSiSj (2-1) 
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where a] and b] are boson creation operator*^ at site i, 
pseudo-spin operator Si = \B\aBi with Bi = {ai^biY, 
a are the Pauh spin matrices, and denotes nearest- 
neighbor (NN) sites of the lattice. We shah consider both 
the square and triangular lattices in the following study. 
Physical Hilbert space of the system consists of states 
with total particle number at each site less than unity 
(the local constraint). In order to incorporate the local 
constraint faithfully, we use the following slave-particle 
representation 



2,3 



i2i 



^i2^i2 - 1 Iphys) = 0, 



(2.2) 
(2.3) 



where (j)i is a boson operator that annihilates hole at site 
i, whereas cpu and (p2i are bosons that represent the 
pseudo-spin degrees of freedom. |phys) is the physical 
state of the slave-particle Hilbert space. 

The previous numerical study of the B-t-J mode]^^^ 
show that there appear various phases including super- 
fluid with Bose condensation, state with the pseudo-spin 
long-range order, etc. For the most of them, the MC sim- 
ulations show that density fluctuation at each lattice site 
is not large even in the spatially inhomogeneous states 
like a phase-separated state. From this observation, we 
expect that there appears the following term effectively. 



Vo 



EH 



Plif + {v\i^2i - P2if 



(2.4) 



where pu etc are the parameter that controls the den- 
sities of a- atom and 6- atom at site z, and Vq{> {)) con- 
trols their fluctuations around the mean values. It should 
be remarked here that the expectation value of the par- 
ticle numbers in the physical state |phys) are given as 

j^HM^i) = ^EiTrphys(4ai) = ^EiTrphys(v^L^ii) 
and similarly j^Y^M^i) = ^ Ei Trphys(^L^2i), where 
N is the number of sites of the lattice and Trphys denote 
the trace over the states satisfying the local constraint, 
i.e., the physical-state condition ()2.3p . Therefore the con- 
straint (j2.3p requires Ea=i /^cri = 1 at each site i. The 
values of Vq and paiicr = 1, 2, 3) are to be determined in 
principle by t, J^, J and filling factor, but here we add 
Hy to Htj by hand and regard parameters in Hy as a 
free parameter. In this sense, we are considering an ex- 
tended bosonic t-J model. It should be stressed here that 
the strong on-site repulsions between atoms enhance the 
effects of spatial lattice and as a result quasi-excitations 
with Lorentz-invariant dispersion can appear in a non- 
relativistic original model like the Bose-Hubbard model^. 

The existence of Hy in Eq. (j2.4p is very useful for 
study of the quantum many-particle systems. In 
the path-integral representation of the partition func- 
tion Z, the action contains the imaginary terms like 
/ dr^i{r)dr(l)i{r), where 4>i is for (j)] and r is the imagi- 



nary time, i.e.. 



= J[D^D(piD(f2]ex.p - J dr(^(pii{r)dr(fii{r) 



^^2i{r)dr(p2i{r) + 0i(r)Sr^i(r) 

-^Hu^HyJ 



(2.5) 



where Htj is expressed by the slave particles and the 
above path integral is calculated under the constraint 
(|2.3p . (In this paper we set = 1.) For the existence 
Hy, we separate the path-integral variables (p^s and (j) as 



^li = VPi7T7iiexp(zcjii), 
= V P^i +-^22 ex.p{iuj2i)j 
(t>i = V Psi +4i ex.p{iujsi), 



(2.6) 



and we integrate out the (fluctuation of) radial degrees of 
freedom. There exists a constraint like in + i2i + ^3i = 
on performing the path-integral over the radial degrees 
of freedom, i.e., i^i {o' = 1,2,3). This constraint can 
be readily incorporated by using a Lagrange multiplier 
A.(r), 



l[S{£u^£2^^i3^) = J 



The variables £cri (cr = 1, 2, 3) also appear in i^tj, but we 
ignore them by simply replacing ip^i y^p^exp(za;cri), 
and then we have 



dXidlie^"^^ EJ=l(-V^o(C^)'+^C^(a.u;^,,+AO) 

dA,e-^^^"^-^^-"--+^^^', (2.7) 



where we have ignored the terms like J drdrOOa^i by the 
periodic boundary condition for the imaginary time. The 
resultant quantity on the RHS of (j2.7p is positive definite, 
and therefore the numerical study by the MC simulation 
can be done without any difficulty. It should be remarked 
that the Lagrange multiplier in Eq. (|2.7p behaves as a 
gauge field, i.e., the RHS of ()2.7p is invariant under the 



following "gauge transformation" , uja^i oJo 



-Oii, Xi 



Xi — drOLi- In the practical calculation, we shall show 
that all physical quantities are invariant under the above 
gauge transformation. 

Various phases can form in the system ()2.ip at T = 0, 
e.g., states with long-range spin orders (FM, AF, spi- 
ral, etc), superfluid with Bose condensation of a and/or 
h atoms, and superposition of them (the supersolid 
(SS))ii^. In the following sections, we shall consider sev- 
eral specific cases of the extended B-t-J model i^t j + Hy 
and derive effective field-theory model. Then we clarify 
the phase diagram of the B-t-J model by studying the 
effective model. 
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III. PHASE DIAGRAM OF THE EXTENDED 
QUANTUM XY MODEL: NUMERICAL STUDY 

A. FM coupling on square lattice 

In this section, we shall study the effective field theory 
obtained in the previous section on the two-dimensional 
(2D) square lattice as the first example. One of the sim- 
plest case is xy-FM system corresponding to = 
and J < in Eq. (j2.ip . though an extension to the 
case \Jz\ < —J is rather straightforward. In this case 
the pseudo-spin symmetry is 0(2), and then we put 
Pii = p2i = P and psi = 1 — 2p. The partition func- 
tion ZqXY of the effective theory is obtained as follows 
from Eqs.(|23|) and (12^]) . 



^qXY 



where 



cr=l,2,3 



i + ^i)^7 



(3.1) 



- ^ Ch COs((^2,i - ^2,j) + COs(f23,i - ^Sj)^ 



(3.2) 



where 
and 



Ch ^p(l - 2p), Cs 4J/. 



(3.3) 



It is obvious that the integrand in Eq. (j3.ip is positive 
definite, and therefore the numerical calculation of ZqXY 
by means of the MC simulation is possible without any 
difficulty. The model (|3.ip is a kind of lattice rotor model. 

For practical calculation, we introduce a lattice for the 
imaginary-time direction and use the following lattice ac- 
tion Alt corresponding to 

3 

^Lr = ^ ^ COs{uJ^^r+r - UJ^^r + K) , (3.4) 
r (7=1 

where r denotes site of the space- time cubic lattice, Cr = 
and At is the lattice spacing of the imaginary-time 
direction. 

We numerically studied the lattice model defined by 
the lattice action 
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FIG. 1. (Color online) Phase diagram of the effective field the- 
ory ^Lxy on 3D cubic lattice for various values of c^. p = 0.35. 
Dots denote the observed phase transition point. Thick line 
denotes first-order phase transition line whereas the others 
are second-order transition lines. System size L — 24. 



where Ai^^e^^"" *^'^) corresponds to A(e*^'^,e *^'^) in 
Eq.([321), 



Ai^{e'^%e-'^-) = 

Cs COs(r^2,r — ^2,rO + COs(r^3^r — ^Sy ] 



-Ci ^ cos(l]i,r — ^ly)^ 



(3.6) 



where Ci = CgAr and Cs = ChAr, and (r, r') denotes the 
NN sites in the 2D spatial lattice. We fixed the value of 
Cr and calculated the partition function Z, the "internal 
energy" E and "specific heat" C as a function of Ci and 
Cs, 



E = (ALxy)/i', 

C = {{A^^y-Ef)/L\ 



(3.7) 



(3.5) 



where L is the linear size of the 3D cubic lattice. In 
order to identify various phases, we also calculated the 
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following pseudo-spin and boson correlation functions, 

i 
i 

Gb{r) = ^(e'"^..e-^"^-*+-), (3.8) 

i 

where sites i and i -\- r are located in the same spatial 
2D lattice, i.e., the equal-time correlations. For exam- 
ple, Ga{r) finite as r ^ oo indicates Bose-Einstein 
condensation (BEC) of the a-atom. 

For numerical simulations, we employ the standard 
Monte-Carlo Metropolis algorithm with local update^^. 
The typical sweeps for measurement is (30000 ^ 50000) x 
(10 samples), and the acceptance ratio is 40% ~ 50%. 
Errors are estimated from 10 samples with the jackknife 
methods. 

We show the obtained phase diagram for p = 0.35 in 
FiglH There are three phases and they are separated 
with each other by first or second-order phase transi- 
tion lines. It is obvious that for small Ci and C3, there 
exists a "paramagnetic state" (PM state) without any 
long-range orders (LRO's) and its domain is decreased 
for increasing c^-. In this PM state, a superposition of 
a-atom and 6-atom is realized at each site, but coherence 
of the relative phase in the superposition does not exists. 
It should be notice that the atomic BEC always accom- 
panies the pseudo-spin LRO. The obtained results are in 
good agreement with the result of the previous study on 
some related model on 3D cubic lattice at finite T-. 



B. AF coupling on triangular lattice 

In this subsection, we shall study the B-t-J model on 
the 2D triangular lattice. As we also discretize the imag- 
inary time, the model ^Lxy in Eq. (|3.5p is defined on the 
3D stacked triangular lattice. In this section we consider 
the case of the xy-AF case with 0(2) symmetry. There- 
fore we set J > and = in Eq. (|2.ip . As there exists 
the frustration, we expect that various phases appear in 
contrast to the xy-FM case. In later section, we also show 
the result of the numerical study of the B-t-J model on 
the stacked triangular lattice at finite but low T, which is 
closely related to the 2D model at T = 0. This close rela- 
tion between 3D model at low T and 2D model at T = 
has been previously observed in various systems^ ^. Ori- 
gin of this similarity is somewhat obvious from in 
Eq. ()3.4p as it can be regarded as an inter-layer coupling. 

To perform the numerical calculation, we first assign 
values of the parameters pu^ p2i {psi = ^ — Pu — p2i)- 
Here we assume a homogeneous state and put pu = p2i = 
p {psi = 1 — 2p) as in the previous case. Results of more 
general cases will be reported in a future publication. 

In Fig|2l we show the obtained phase diagram for 
p = 0.3 and the hole density= 0.4. There are six phases 
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FIG. 2. (Color online) Phase diagram of the quantum XY 
model on triangular lattice. Density of each boson p = 0.3. 
System size L = 18. There are six phases. Phase of 120^ de- 
notes the state with 120° pseudo-spin long-range order with- 
out atomic BEC. SF stands for superfluid. 

in the Cs — Ci plane. For small C3, the pseudo-spin 
degrees of freedom exhibits the 120^ long-range order. 
As we assume a homogeneous hole distribution, a super- 
posed state of atoms and hole is realized at each site, but 
coherent condensation of atoms does not takes place yet. 
As C3 is increased, phase transition to the states with 
the 120^ long-range order and Bose condensation take 
place. For example in the phase B in FiglSl the both 
atoms Bose condense, as the correlation functions shown 
in Figl3] exhibit. In the mean-field approximation, the 
wave function of that state is given by 

\'i^)=Ui4 + bl + CA] 

X Ylia] ^e'^bl ^cb] 
ieB 

iec 

where c's are some complex number and |0) is the empty 
state of a and 6-atoms. This phase will be discussed by 
using the effective field theory in Sec. IV. C and also by 
the MC simulations in Sec.V. As the parameter C3 is 
increased further, order of the pseudo-spin is destroyed 
first (phase C) and then changes to the FM one in the 
S"^ -Sy plane (phase D). 

IV. EFFECTIVE FIELD THEORY FOR BOSE 
CONDENSATION AND PSEUDO-SPIN ORDER 

In the previous sections, we numerically studied lat- 
tice quantum XY model that describe the two-component 
bosons with strong repulsions. The obtained phase dia- 
grams show that there exist Bose-condensed phases as 
well as the state of simple pseudo-spin LRO. In this sec- 
tion, we shall derive an effective field theory from the 
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and the fact that other correlators Hke (^^^^ii^) ^^^ii^ )^ 



Gs(r) _ 

Ga(r): 

Gb(r): 



FIG. 3. (Color online) Correlation functions in the phase 
diagram. C3 = 0.65 and Ci = 1.6 (left), C3 = 1.0 and Ci = 
1.6 (right). 



quantum XY model, which describes directly the Bose 
condensation and spin order. This field theory not only 
explains the numerically obtained phase diagram but also 
reveals low-energy excitations and interactions between 
them. This kind of effective field theory has been ob- 
tained and discussed for a granular superfiuid etc by 
using Hubbard- Stratonovichi transformation^^. In this 
section we shall employ a slightly different method^. As 
a result, we reveal some important point that has been 
overlooked so far. 



A. Effective field theory 

To illustrate the procedure to derive the effective field 
theory, we first consider a single rotor model as an ex- 
ample that describes superfiuid phase transition. After 
consideration of the simple model, we shall apply similar 
methods to the quantum XY model for the B-t-J model. 

We first rewrite the partition function of the rotor 
model on the square lattice by introducing source terms, 



=/[. 



\uj\e 

[dUj\e ^ X J ^Srji^ SrjiJ 

Xl{r]i,UJi)\r^=f^=0, 
/(7/„C^,)=e^'"^^"^'^"+^^'"^"\ 

where = drUJi and 

(id) 

In Ea. (j4.ip . we evaluate the path integral of as 



(4.1) 
(4.2) 



g/dr/dr'f7.(r)e-^ok-'l^^(^')^ ^^^^^ 



where we have used 



are vanishing. RHS of Eq. (|4.3p can be expressed by in- 
troducing a complex boson field as 



1 



J [d^] exp 



(4.5) 



By inserting Eq. (j4.5p into Eq. (j4.ip . we obtain the effec- 
tive field theory of the rotor model with the following 

action Arotor, 



dr 



^2 + 



-^^($:$,+c.c.) 



(4.6) 



The above derivation from Eq. (j4.ip to Eq. ()4.6p seems 
exact because the integration of 0Ji{r) is the one-site inte- 
gral and essentially Gaussian integration of the free fields. 
In fact the above manipulation is exact as long as the bo- 
son field does not Bose condense. On the other hand 
for J Vo^ the Bose condensation of takes place, i.e., 
{^i) 7^ 0. One may wonder if the above manipulation 
is applicable even for this case because Eq. ()4.4p seems 
to indicate nonexistence of the long-range oder of 
Furthermore in this case, the mass term of becomes 
negative, and the integration of in Eq. (|4.6p becomes 
unstable, e.g., for the square lattice 

J 



-^(^*^,+^ + c.c.)-FoEl^^l 



(4.4) 



2,/X I 

= -^^Y.\'^,H'' + Y.i^j-v,)m\ (4.7) 

where is the lattice difference operator and z is the 
number of links emanating from a single site and z = 4 for 
the square lattice. This instability comes from the fact 
that the order of the 77-derivative and the cj-integration 
is not interchangeable when the Bose condensation, i.e., 
a phase transition to an ordered state, takes place. It is 
obvious {^i) = (e*'^*) < 1 in the Bose condensed state, 
and therefore it is plausible to expect term like A|$i|^ to 
appear in the effective field theory to stabilize the inte- 
gration of for the Bose condensed state, though ex- 
plicit calculation to determine the coefficient is difficult. 
Simple spin- wave like approximation for the rotor model, 

JcOs(cc;^ — UJj) ^ J — ^{^i ~ 

Vq 

gives the estimation like (e*^^)(e~*^^) ~ e~~, and then 
the coefficient of the |<l>^|^-term is estimated as A ~ Je~ 
for J> Vb. 
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Hubbard- Stratonovich transformation derives a similar 
effective field theory to the above. But its straightfor- 
ward application gives a negative coefficient of the 
term indicating an instability of the system. This means 
that certain step of the derivation is invalid, e.g., intro- 
duction of the Hubbard- Stratonovichi field and the uj- 
integration is not interchangeable. This problem is un- 
der study and the result will be reported in a separate 
publication. 

Similar manipulation to the above can be applied 
to the present quantum extended XY model. From 
Eqs.dSD) and (jT^ . 



V = 



^qXY = / [duj]e 



[d(jj\e ^ az ^Sv^i^Sv^i^ 

^ (j](Ti') ^ai) 177=77=07 



(4.8) 



where u^ji = Sr^ai and we have omitted the gauge field 
as we consider the only gauge- invariant objects through 
flat- Then the integration over can be performed as. 



I{r]ai,^ai) = HVai, Vai), (4.9) 



and the Green function of ujai{r)^ G^j^ is obtained as 

^eic^^.(r)g-iu;./.(r')^ = 5^^/ 6"^° l^-^'l . (4.10) 

Typical term of I{r]ai^f]ai) is as follows, 

and this quantity is expressed by introducing complex 
scalar fields ^aii'T') as 



J [d^] exp 



1 



/ 



4Fo 



(4.11) 



Ai|$i|4 + 2A2|$2|^ 

+A3|$lP + 2A4|$2|^ - ^|$l||$2|^ 




01 = 1.02771 02=1.11237 
V= -2.05918 



FIG. 4. (Color online) Potential V{^a) (|4.14p and its mini- 
mum. Symmetric case $2 = $3. Parameters are Ai = 1, A2 = 
^,A3 = — 1,A4 = —1, and g — 1. 



system. The final form of the effective action is therefore 
given by 

-I' 

Z = J [d^]e 



eff, 



(4.13) 



By inserting I{r](ji,f]cri), which is expressed in terms of the 
boson fields ^aii^) for e*^'^*, into Eq. (j4.8p . the action of 
the effective field theory is obtained as follows, 

Ao = J dr[ ^ (a.$;^$,^.)-^E(l*-l' 



B. Phase diagram and low-energy excitations: 
Square lattice 

In this subsection, we shall apply the field-theoretical 
approach explained in the previous subsection to the B- 
t-J model on the square lattice as a simple example. Fur- 
thermore we focus on the FM parameter region J < 0. 
In this case there exists no frustration and therefore it is 
rather straightforward to obtain the phase diagram. Nev- 
ertheless study on this system reveals important aspect 
of the state with multiple long-range orders and structure 
of the Nambu-Goldstone bosons. 

This system was studied in Sec. HI. A, and we obtained 
the phase diagram by MC simulations. The potential 
V{^a) of the present system is given as follows from the 
effective field theory in Sec. IV. A, 

V{^a) = ^{Vo - 2za^)\^^\^ - g{^l^2n + c.c.) 



-E5($ti*2i*3i + C.C.) 



(4.12) 



where ai = a^, a2 = as = and g = To describe 

the Bose condensed state, we add the |<I>|^-terms in the 
effective field theory and discuss the phase diagram of the 



From the first term of Eq. (j4.14p . it is obvious that spon- 
taneous symmetry breaking occurs for > Vq, but the 
second terms that represent the interplay between the or- 
der parameters give nontrivial contribution to the phase 
diagram. In FigHl we show typical potential and its 
minimum obtained from Eq. (|4.14p . which derives qual- 
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itatively the same phase structure with that obtained in 
Sec.IILA. 

It is interesting to see how gapless low-energy excita- 
tions, i.e., Nambu-Goldstone (NG) bosons, appear in the 
present system. In the phase with the FM spin order 
{^s) = V > 0, we set = v -\- i/ji -\- ixi^ and then it is 
obvious that the field Xi describes a NG boson that ap- 
pears as a result of the spontaneous symmetry breaking 
of the U(l) pseudo-spin symmetry. 

Interesting point is that how many NG bosons appear 
in the phase with the FM+2SF. If ^ = in Eq. (|4:T2|) . it is 
obvious that there exist three NG bosons. However in the 
original B-t-J model, the symmetry is U(l)xU(l) for the 
global phase rotation of a and 6-boson operators. In order 
to study the low-energy excitations in the effective field 
theory, we shall take the continuum description instead 
of the lattice one though it is not essential. 

For the FM B-t-J model on the square lattice, it is 
straightforward to derive effective Lagrangian Leff in the 
continuum spacetime from the effective field theory on 
the lattice. For example in Eq. (|4.12p . we put 



where 



Then Leff in the continuum is derived as follows from 

Eo. Km 



^eff 



1 



-^|a.$3(a:)n+n*<.), 



(4.15) 



where the potential V(i>cr) is given as 

= {Vo - 2za,)|$i(x)|2 + {Vo - 2zah)\Mx)f 
+iVo-2zah)\<i>3ix)f 

-g(^l{x)^2ixMix)+c.c.) (4.16) 
+Ai|$i(x)|4 + X2\Mx)\^ + A3|*3(a;)|^ 
We first substitute 

into V{^^), 

V(no, ni) = Aino + 2X2nl + X^uq + 2X^711 

-2^V^ni, (4.17) 

and derive the equations that determine the value of 
no, ni as, 

^ 2Aino + A4 - 0, 

dV 
dni 



= 4A2ni + 2A5 - 2g^ = 0, (4.18) 



A4 = Vb - 2zas, 
A5 = Vb - 2zah. 



(4.19) 
(4.20) 



Next, we introduce three phase fiuctuation fields that 
represent NG bosons in the FM+2FS phase with no 7^ 
0, ni ^ 0, 



^i{x) = + i(/)(x), 
^2{x) = ^/n^ ^ ibi{x), 
^six) = ^/rH ^ ib2{x). 



(4.21) 



By substituting Eg. dOT]) into V(<^>a) in Eq.lQJ^, the 
quadratic terms of the fluctuating fields are obtained as 

2Aino + A4 -g^/n{ g./rH 
^2 = (0,^1,^2) I -gVrii 2A2ni + A5 -g^/rl^ 
gJW -g\/no 2A2ni + A5 



b2 



By using Eq. (|4rT8l) . 



V^2 = ((/>, 61, 62) 



/no -dV^ 

-gV^ gV^ -gV^o] I 

-gx/n^ gJno I \ b2 



hMM)g^ 




(4.22) 



It is straightforward to diagonalize the matrix K by using 
a unitary matrix U and obtain the mass gap. 



^0 

U-^KU= I 
.0 




(4.23) 



Therefore there are two gapless modes that correspond to 
the NG bosons, though one may expect three NG bosons 
as the U(l) pseudo-spin symmetry is spontaneously bro- 
ken and also both the a and 6-atoms Bose condense. 
From the above derivation of the mass gap, it is obvi- 
ous that the cubic coupling g^l^2^s plays an essentially 
important role. 



C. Phase diagram and low-energy excitations: 
Triangular lattice 

In this subsection we shall focus on the system on the 
triangular lattice and study the phase structure and low- 
energy excitations. The obtained results will be com- 
pared with those by the numerical study on the quan- 
tum XY model given in the previous section and those of 
the 3D B-t-J model in the following section. To study 
the system, we assume a homogeneous state and put 
Pii = P2i = P {p3i = 1 - Pli - P2i)' The field $1 repre- 
sents pseudo-spin degrees of freedom, and it is plausible 
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to assume that its condensation has a uniform amphtude 
=constant as holes are distributed homogeneously. 
On the other hand, its phase degrees of freedom has a 
nontrivial behavior as a result of the frustration. More 
general cases will be studied in a future publication. 

We first study the effective action by a mean-field the- 
ory like approximation assuming the a/3 x ^/S symmetry. 
The triangular lattice is divided into three sublattice la- 
beled A, B and C sublattices. Each field on the sublattice 
A is denoted as etc. It is not so difficult to search 
the ground state of the potential. For small ah (i.e., small 
t), <l>2 and $3 do not Bose condensate. On the other 
hand, a nontrivial long-range oder of the pseudo-spin ap- 
pears for tts > Vq. For constant |<l>i|, one can show that 
the state with the three sublattice coplanar order like 
^lA = l^il, ^iB = |^i|e^2V3^ = |$i|e-^2V3 ap- 

pears as the lowest-energy state. This state obviously 
corresponds to the phase of the 120^ spin order in Figl2l 

To obtain the expectation value of and to study 
low-energy excitations, we relabel the lattice site by 
dividing it into three sublattices as z (^5, o) where 
o = A^B and C. Then we parameterize the field $i 
as 



^ia{s) =pi 



(4.24) 



Effective Hamiltonian of the spin part H^^^^ is readily 
obtained from Leff in ()4.13p (please notice < in the 
AF case), 

^f:r = E(i«^i*H*i^+^+^-c 

+ E(^o|*ii|' + A3|$iin, (4.25) 
and substituting Eq. (|4.24p into i^^ff 

^eff" = ^COS^Mas|piPi+^ +C.C 



+ E(^il'°*l'-(3k|-^o)|Pin,(4.26) 

where denotes the difference operator on the trian- 
gular lattice. 

From Eq. (j4.26p . Bose condensation of pi takes place 
for 3as > Vo, and the classical expectation value pd of pi 
is easily obtained as 



Pel 



'3\as\-Vo 



2Ai 



(4.27) 





FIG. 5. (Color online) Dispersion relation ^(k) in Eq. (|4.33p . 

Parameters are as = l,Vb = l,7 = 5 and {3—1. 



To study the low-energy excitations, we introduce a 
complex scalar field rji as follows. 



Pi = Pel + ^i, 



(4.28) 



and then from Eq. (j4.25p . we obtain the quadratic terms 
of rji in the Hamiltonian as follows 



(4.29) 



where 



X(k) 



-2|a.sl cos I — I cos I — - — 



— \as\ cos kia + 7, 
f3 = Xip% 7 = 4Aip^^ + K). 



(4.30) 
(4.31) 



We diagonalize Eq. (|4.29p by the Bogoliubov transforma- 
tion and obtain 



where 



^(k) = Vx(k)' - 4/32, 



(4.32) 



(4.33) 



and its typical behavior is shown in FigEl As in the limit 



E{k) 0, 



(4.34) 



then 6k represents a NG boson corresponding to the spon- 
taneous symmetry breaking of the U(l) pseudo-spin sym- 
metry. 
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Finally let us study the low-energy excitations in the 
phase with the 120^ spin order plus 2SF's that corre- 
sponds to the phase B in Fig 121 From Eq. ()4.13p . the ef- 
fective Lagrangian L*^ of low-energy excitations in this 
phase is obtained as 



(4.35) 

^(^-) = E (^0 - 2z\as\^ (Vo - 2zah^ |$2,^|' 



Vo-2zah]\^: 



Ail^i 



(4.36) 



A2l<^>: 



3,i 



where z is again the number of links emanating from a 
single site. 

To obtain the ground state, we adopt the following 
Ansatz for the various condensations, which respects the 
three-sublattice symmetry. 



^2,A Vth, ^2,b = Vrne'^, ^2,c = yjn{e~'^ , 



3, A 



ni,<l>3,. 



,^3,c = Vme'^. (4.37) 

The values of no, n\ and /3 are determined by substitut- 
ing Eq. (j4.37p into (j4.35p and imposing stationary condi- 
tion, 

Ql}^'^ /27r \ 

UahTii sin/3 - 4.g^ni sinf y - 2/3 j = 0, 



^^2Aino + A4-^(5^ 



dni 



0, 



4A2ni + 2A5 - 2gSy/n^ = 0, 



(4.38) 



where (5 = |(1 + 2cos(^ - 2/3)). It is verified that non- 
trivial solutions of no and ni exist for sufficiently large ah 
and as (i.e., negative A4 and A5). Correlation functions 
obtained from a typical solution to Eq. ()4.38p are show in 
Figini which have similar behavior to those obtained by 
the previous MC simulations of the quantum XY model. 

As in the square lattice case, we introduce fields that 
describe low-energy excitations. 



^l,A = 


V 










{^ 


/no- 


f i(/))e*^. 




^1,C = 


(. 


/no- 


f i0)e"^^. 


(4.39) 


^2,A = 


V 


^4 


- ihi, 




^2,B = 


(' 


s/rii 






^2,C = 


(^ 


\/rii 




(4.40) 



SF con-elations C3=1.0 Cl=1.5 



Effective ■ 
MC : • 



FIG. 6. (Color online) Boson correlation functions G^{r) — 
($2,i^2,i+r) = ($3^^$3,i+r) in the effective field theory, and 
Ga{r) = Gb{r), which is numerically obtained in the quantum 
XY model. We set no = ni = 0.5, g — 6, as = 1 and ah = 1. 



^3,A = Vni^ib2, 
^3,3 = (v^ + ^^2)e"'^, 



(4.41) 



Substituting the above equations (|439l) . (|0Q1) and (|OT]) 
into L*^ in ()4.35p . we obtain the mass term of the fiuc- 
tuating fields as 



g^-7^ 

V2 = (0,61,62) -gSJnl g5Jn^ 



g8Jn{ 



gS^/TH g^^/nl 

g^Vno 
g5y/n^ g5y/n^ 



{^,bi,b2)gSK^ I bi 
b2 



(4.42) 



where we have used Eq. (|4.38p . We can easily diagonalize 
and get eigenvalues. 






Q Q 2no+ni 



(4.43) 



The above result indicates that there are only two NG 
bosons, though the pseudo-spin U(l) symmetry is spon- 
taneously broken and both the a and 6-atoms Bose con- 
dense. 



V. PHASE DIAGRAM OF THE BOSONIC t-J 
MODEL ON STACKED TRIANGULAR LATTICE: 
NUMERICAL STUDY 

In Sec. III. B and IV. C, we studied the low-energy effec- 
tive theories of the B-t-J model on the triangular lattice 
at T = and obtained phase diagram and low-energy 
excitations. As we argued previously, if fluctuations of 
particle density at each site is not so large, the model at 
T = reduces to the system with "Lorentz invariance" , 
i.e., a linear-time derivative term change to a quadratic- 
time derivative term^. In this case, the imaginary time 
r plays a similar role to the inter-layer dimension of the 
stacked 3D lattice. 

In the present section, we shall study the B-t-J model 
on the stacked triangular lattice by numerical simulations 
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FIG. 7. Phase diagram of the bosonic t-J model in a stacked 
triangular lattice, pa = = 0.3, and ci = 10.0. 
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FIG. 9. (Color online) Spin and boson correlation functions 
obtained by MC simulation of the B-t-J model in a stacked 
triangular lattice, — 15.0 and C3 = 5.0 in Fig|7l Css = 

where sites i and i -\- r are located in the same 2D triangular 
lattice. The result indicates the phase with 120"^ spin order 
and 2SF's. 

previous studies in this paper, we employ the canonical 
ensemble with the particle number of each atom fixed. To 
impose the local constraint, we employ the slave-particle 
representation (j2.2p . Then the partition function Z is 
given as 




FIG. 8. (Color online) Snapshots of pseudo-spin (upper 
panel), density of a-atom and 6-atom (lower panel). — 
is the easy plane of the pseudo-spin (see Fig|9|), and length of 
arrows indicates magnitude of the pseudo-spin. A stripe oder 
forms and also there are hole-rich regions. 

with a quasi-classical approximation. Hamiltonian on the 
stacked triangular lattice is given as follows, 

^3DtJ = - ^ t{alar' + blbr' + h.c.) 

(r,r') 

- t\ala^^^ + blb^^^ + h.c.) 

r 

+ jY,iSrSr'+S?Sf,) (5.1) 

(r,r') 

where (r, r') denotes the NN site of the 2D triangular lat- 
tice, and 3 is the unit vector in the inter-layer direction. 
There are (at least) two ways of the MC simulation, i.e., 
the grand-canonical and canonical ensemble. As in the 



Z = J [D(t)DipiDif2]e-^^^''''. (5.2) 

The path integral in (|5.2p is performed by the MC sim- 
ulation with local update keeping each particle number 
fixed. We call the calculation (|5.2p the quasi-classical ap- 
proximation as we ignore the Berry phase in the action 
of the path integral. Some detailed discussion on the va- 
lidity and physical meanings of this approximation was 
given in the previous papers 

We first show the phase diagram for = P6 = 0.3 and 
the hole density = 0.4 in Figlll where ^ = t' /{k^T), c'^ = 
J / {k^T) and C3 = t/{k^T). For small hopping amplitude 
t, the system forms a solid with voids whose snapshot is 
shown in FiglHl As C3 is increased, phase transition to 
a state with the 120^ spin order and 2SF's takes place. 
In Figini we show the spin and particle correlation func- 
tions that verify this conclusion. The previous numerical 
study of the quantum XY model and the analytical study 
by the effective field theory predict the existence of this 
phase. This result again indicates a strong resemblance 
of the phase diagram of 2D system at T = and that of 
the corresponding model in stacked 3D lattice. As C3 is 
increased further, phase transition to a phase-separated 
state takes place. In this phase, the system is divided 
into a-atom rich region and 6-atom rich region, and in 
each region a SF of single atom forms. 
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VI. CONCLUSION 

In this paper we studied the B-t-J model of the two- 
component bosons with strong on-site repulsions. We 
used the salve-particle representation to treat the local 
constraint faithfully and derived the effective field theory 
by integrating out radial degrees of freedom of the slave 
bosons. The resultant field theory is a kind of extended 
quantum XY model, and its phase diagram was investi- 
gated by the MC simulations. The results shows that the 
B-t-J model on the square lattice has a simple phase dia- 
gram whereas that with the AF pseudo-spin coupling on 
the triangular lattice has rather complicated structure. 

Then we derived the second version of the effective 
field theory of the B-t-J model by using a "Hubbard- 
Stratonovich transformation" with the source terms. The 
resultant field theory describes the pseudo-spin and Bose 



condensation of the a and 6-atoms directly. The effective 
potential of these order parameters clarifies the phase 
diagram of the system and verified the existence of the 
interesting phases like 120^ (pseudo-spin) +2SF's on the 
triangular lattice. Furthermore, we studied low-energy 
excitations, in particular the NG bosons, in the effec- 
tive field theory and found that there exist only two NG 
bosons even in the phase in which the U(l) symmetry of 
the pseudo-spin is spontaneously broken and the both a 
and 6-atoms Bose condense. 

Finally we investigated the B-t-J model on the stacked 
triangular lattice at finite T by the MC simulations and 
showed that a similar phase diagram to that of 2D system 
at T = appears. This result means that the direction 
perpendicular to the 2D triangular lattice plays a similar 
role to the imaginary-time direction because of the strong 
on-site repulsion and the spatial lattice structure. 
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